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Cavitation Density of Superfluid Helium-4 around 1 K 
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Using an optical interferomertric method, the homogeneous cavitation density of superfluid helium 
at T = 0.96K is measured and found to be p cav = 0.1338 ± 0.0002g.cnU 3 . A well established 
equation of state for liquid helium at negative pressures converts this to the cavitation pressure 
P cav = —5.1 ± 0.1 bar. This cavitation pressure is consistent with a model taking into account the 
presence of quantized vortices, but disagrees with previously published experimental values of P ca v 


I. INTRODUCTION 

Because helium is a model material at low tempera¬ 
ture, its phase transitions have been studied in much de¬ 
tail. In particular, the stability limits of the liquid phase 
with respect to the solid phase at high pressure and to 
the gas phase at negative pressure have been studied ex¬ 
perimentally and theoretically (see the review articles of 
S.Balibar et aZ.p],0). To avoid heterogeneous nucleation 
of the new phase, experimental over- and under-pressures 
are produced in bulk liquid away from any surface using 
high amplitude focused sound waves. While the appear¬ 
ance of the new phase is easily detected optically, the 
measurement in situ of the local pressure at which nu¬ 
cleation occurs is a challenge. Estimations have been 
drawn either from the oscillation amplitude of the sound 
emitter combined with an estimation of the gain due to 
the focusing , or from extrapolating the liquid static 
pressure down to where a vanishingly small sound wave 
would produce cavitation [h]. However nonlinear effects 
in the sound wave propagation introduce uncertainties in 
both methods. 

In 2010, our group introduced a time-resolved quan¬ 
titative multiphase interferometric imaging technique [Q] 
for measuring the density of a medium inside a sound 
wave with cylindrical symmetry. In 2012, this tech¬ 
nique was successfully applied for density measurements 
in solid helium at 1 K in pressure swings below the equi¬ 
librium melting pressure Q ■ We intended to implement 
the same technique to study metastable liquid at pres¬ 
sure above this melting pressure. In order to check the 
method in liquid helium, we decided to measure first the 
liquid density at which cavitation occurs around 1 K aim¬ 
ing at verifying the results of reference j5|. In this article, 
we present our result for the cavitation density. Previous 
experimental results were given as cavitation pressures. 
So an equation of state for liquid helium at negative pres¬ 
sures have to be used for comparison. Fortunately, var¬ 
ious theoretical approaches |8Hllj have produced quite 
similar equations of state. To our surprise, our density re¬ 
sult converted to pressure does not quite agree with pre¬ 
vious estimations. Thus we also reproduced the pressure 
extrapolation of F.Caupin and S.Balibar Q, and found 
a reasonable agreement with their data. After discussion 
of various sources of uncertainties in our measurements, 
we compare our results with various theoretical estimates 


of the cavitation pressure. 

II. EXPERIMENTAL SET-UP 

The technique to measure the local density variations 
in focused acoustic waves has been described in previous 
articles @, [7]. Here is a brief reminder. 

The experimental cell containing liquid helium is 
cooled in the cryostat with four optical ports. The work¬ 
ing temperature can be regulated from 0.9 K to 2.1 K. In 
the experiment described in this paper, the temperature 
is fixed at T = 0.96 K. The cell is connected to a buffer 
volume at room temperature, so that the static pressure 
P s t is easily monitored with the help of a Keller X35 pres¬ 
sure sensor with an accuracy of ±15 rnbar. A hemispher¬ 
ical piezoelectric transducer (PZT) excites and focuses 
ultrasound waves in helium at the frequency 1.15 MHz of 
its first thickness vibration mode. The transducer inner 
diameter is 12 mm, and the thickness is 2 mm. In order to 
observe directly the acoustic focus, a small part, 0.9nnn 
in height, has been removed around the transducer rim. 
One side of the transducer is grounded and the other side 
is connected to the output of a RF amplifier driven by an 
arbitrary function generator (AFG). A detailed scheme 
is shown in FigjT] 



FIG. 1. Experimental set-up. At the center of the transducer, 
the small star is the acoustic focus where cavitation occurs. 
CCD camera monitors the focus plane through a lens. The 
camera image is composed of pixels distributed in lines and 
columns. We optimized the camera orientation so that lines 
are parallel to the a>axis. 

The optical part of the set-up is a Jamin interfer- 
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ometric imaging system with a pulsed Nd:YAG laser 
(A 0 = 532 nm) as a light source and a CCD camera for de¬ 
tection. The focal region is imaged onto the CCD camera 
with a magnification factor of 2 using an appropriate lens. 
The spatial resolution of the entire system is 20 /jm [f|, 
about 1/10 of sound wavelength A s = 0.24 mm. One arm 
of the interferometer passes through the acoustic focal 
region while the other crosses the cell in an unperturbed 
zone. The acoustic wave introduces a density variation 
that modulates the refractive index of helium (mainly 
in the acoustic focus) and hence gives rise to an optical 
phase shift between the two paths. A pair of phase plates 
is placed between the laser and the cell. One of these 
phase plates is mechanically controlled by a computer in 
order to add a known phase shift to the unperturbed ray. 
This added phase shift enables us to extract the phase 
shift due to the variation of the refractive index as we 
will see in the next section. 


III. LOCAL DENSITY MEASUREMENT 

Using the AFC, we can choose the time t with respect 
to the sound pulse triggering at which an image of the 
interference field is taken by the camera. This time is ad¬ 
justed by steps of 0.05 p, s, about 6% of the sound period. 
By repeating the measurements and recording images at 
successive delays on a relatively long period, we are able 
to reconstruct a temporal evolution of the interference 
field in the focal region. 

The observed intensity of a given pixel at time t is 
expected to be: 

lit) = I 0 (t)[ 1 + C(t) cos (6<j>(t) + /?)] ( 1 ) 

8<f>{t) is the optical phase to be measured, and (3 is the 
controllable phase shift due to the phase plates. Io{t) and 
C{t) are respectively the mean intensity and the fringe 
contrast. A series of measurements with different /3 at the 
same time enables us to extract the phases 5(j>(t) through 
a fit@] of eq.((T]). 

Let y be the light propagation axis, z the PZT axis and 
x the axis orthogonal to y and 2 (FiglTJ . At a given delay 
i, the phase shift S(f>{x , z) is related to the refractive index 
map by a simple integration along the y-axis. In our 
case, the sound field is rotationally invariant around the 
hemisphere axis z, so that the refractive index variation 
5n is only a function of z and r = \J x 1 T y 2 . Given the 
fact that Sn is 0 outside sound field, we can write 5<fi as 
the Abel transform of Sn: 

o_ /*+oo 

8<j>(x, z) = — f dy Sn(\/x 2 +y 2 ,z) ( 2 ) 

Conversely, radial refractive index profiles can be re¬ 
trieved from phase shift maps via an inverse Abel trans¬ 
form. Then, using Clausius-Mossotti relation in the limit 
no ~ 1, the density variation Sp of the medium can be 


deduced easily: 

An_ = ( 3 ) 

n 0 - 1 po 

no and po being the unperturbed refractive index and 
density. 


IV. DETERMINING THE CAVITATION 
VOLTAGE 

Our measurement of local density by interferometry 
requires a completely reproducible phenomenon. If bub¬ 
bles appear randomly, it is impossible to measure the 
optical phase shift 6<f>(t) for the pixels involved. In other 
words, this method only allows us to measure the local 
density just below the cavitation density where no cavita¬ 
tion process occurs (or the cavitation probability is very 
low). Then the measured local density should be very 
close to the real cavitation density and a linear extrap¬ 
olation to the cavitation voltage would introduce only a 
small correction. 

Thus, before performing any density measurement, we 
have to precisely determine the cavitation voltage. Other 
groups have observed that bubble life time in superfluid 
helium depends on the static pressure P st and is of the 
order of some tens of microseconds I,EH. Then, ~ 10 ps 
after the minimum pressure wave front passed the acous¬ 
tic focus, bubbles have expanded to their maximum size 
and are easily observed on the CCD camera (see Fig(2]). 

The cavitation process has a statistical behaviour be¬ 
cause of the thermal fluctuations. According to F.Caupin 
et al. [H , this probability is described by the “asymmetric 
S-curve formula”: 

E(V) = 1 - exp[- In 2 exp(^(I//V r c - 1))] (4) 

where V is the excitation voltage, V c the cavitation volt¬ 
age and £ a dimensionless parameter. In order to deter¬ 
mine V c , we proceed as follows. For a given static pres¬ 
sure, the bubble probability is determined for 5 different 
excitation voltages. Each voltage point corresponds to 
1000 trials (1000 sound pulses) and the probability is then 
simply given by the number of positive events (creation 
of bubble) divided by the number of trials. The relative 
standard deviation on the probability is l/\/T000 ~ 3%. 
To avoid heating, 10 bursts of 100 sound pulses at 10 Hz 
repetition rate were shot, waiting 100 s between each 
burst. To precisely control the driving voltage of the 
PZT, we fixed the RF amplifier gain factor at about 390, 
and adjust only the AFG voltage amplitude with a rela¬ 
tive accuracy of 10 -4 . Hence we use the AFG voltage V 
as a scale to determine the cavitation voltage instead of 
the PZT driving voltage. 

The cavitation voltage is the value corresponding to a 
bubble probability of 1/2 according to eq|U As it can be 
seen in Fig[2l the relative width of the curves is about 
1% of V c . These curves are indeed very sharp and for the 






3 



FIG. 2. Left: Cavitation probability at 0.96 K for three differ¬ 
ent pressures P a t : circles 0.15 bar, triangles 0.65 bar, crosses 
1.26 bar. The corresponding lines are fits according to eq. 0. 
Right: Image of a bubble recorded by the camera 10/xs after 
its creation. 


AFG voltage V max of about 2% below V c , the probability 
Z(V max ) is about 1CT 3 . 


V. CAVITATION DENSITY 


extraction process of the optical phase shift induced by 
the acoustic wave, and the systematic errors arising from 
an imperfect cylindrical symmetry of the pressure wave. 

As we mentioned earlier, for each time and for each 
pixel, the phase shift value 5<j) (t) is obtained by applying 
a fit on /3-dependent intensities using eq. © We use 

a computer program to extract these phases with 95% 
confidence bounds. This gives the phase shift uncertain¬ 
ties mainly due to shotnoise, camera reading noise and 
laser power fluctuations. Once the phase shift map is 
determined, an inverse Abel transform Q is applied to 
recover the refractive index local variation induced by the 
acoustic wave. Then, the density variations are deduced 
from the optical index variations using eq.®. The in¬ 
verse Abel transform is a linear transformation. For a 
given line of an image, the calculation of the optical in¬ 
dex at a given pixel i depends linearly on the phase shift 
values for all pixels on the same line. The local density 
variation at this pixel is thus in the form : 


In this section, we present our results for the cavitation 
density of superfluid helium at T = 0.96 K. The cavita¬ 
tion density was reached from three different static pres¬ 
sures : 0.15 bar, 0.65 bar and 1.26 bar. For each static 
pressure, the cavitation voltage is determined following 
the approach depicted previously (see Fig®. Then the 
minimum density in time and space p m in was measured 
for several voltages below the cavitation threshold. An 
example is shown in Fig® The density variation is not a 
linear function of the voltage. This is due to nonlinear ef¬ 
fects in the acoustic wave and possible to the appearance 
of shock waves jl.3l j. Assuming a local linear dependence 
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FIG. 3. Experimental measurements of minimum densities 
pmin for different driving voltages V at T = 0.96 K and P s t. = 
0.15 bar. The dashed line represents the cavitation voltage. 

of Pmin to V, and taking into account the error bars, one 
can safely consider p cav is PminiVmax)■ In order to im¬ 
prove the accuracy in determining p cav , we have made a 
lot of measurements of the local density just below the 
cavitation voltage. 

The measurement uncertainties can be divided in two 
parts : the statistical ones which come mainly from the 
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5pi ='^2a i jS4>j (5) 

j=i 


where j is the pixel index and a,j is a weight. The er¬ 
rors A Spi on 8pi can be computed from the error A 8(f>j 
on 1 5<f>j and the weights onj which could be in principle 
extracted from the Abel inversion program. Instead we 
used a simpler empiric method, assuming that the phase 
uncertainty is about the same for each pixel, and is not 
correlated from one pixel to an other. In that case, the 
uncertainty of the density variation at the pixel i is : 

OO OO 

(A Sp z f = Y a% (ASfij) 2 = (Ac5</>) 2 Y ( 6 ) 

j=i j=i 

Then, we performed N = 1000 density calculations, for 
the same treated line while adding a gaussian noise to 
the phase shifts for every calculation. The standard de¬ 
viation of the added noise is chosen to be the same as 
the phase shift uncertainty A6<j). Once we have these 
N treatments, the statistical uncertainty (A 5pi) N of the 
radial density variation is calculated for each pixel of the 
line. By construction, this uncertainty is equal to y/2 
times the original unknown statistical uncertainty A Spi 
of the density variation, because : 


(A Spif N = Y‘ 


{AScj)) 2 + (A 5cj)) 2 = 2 (A Spi) 2 


(7) 

Appling this method, we found that the statistical un¬ 
certainty around the cavitation density is on the order of 
0.0002 g.cm -3 , while the maximum value of Sp is of order 
0.0125 g.cm" 3 . 

The inverse Abel transform assumes that the symme¬ 
try axis is exactly known. Actually, it is unknown and 
has to be determined experimentally by searching a sym¬ 
metry axis in the phase maps. But the phase noises as 
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well as any possible asymmetry of the acoustic wave lo¬ 
cally perturb the left-right symmetry of the phase shift 
maps. This perturbation will add an uncertainty in the 
calculation of density variations. The difference between 
the Abel inversion applied to the left and to the right of 
this axis gives an order of magnitude of this uncertainty. 



FIG. 4. Cavitation density as a function of static pressure 
at T = 0.96 K. For more clarity, the different measurements 
are shifted from their actual P„t values (black circles). The 
horizontal continuous line represents the average cavitation 
density and the dashed lines its uncertanty. 

The symmetry axis for a given phase map is found 
by fitting a straight line through all symmetry centers 
when the amplitude of the sound pulse at focus is max¬ 
imum. Then the mean and the standard deviation for 
the position of these axis are computed.The uncertainty 
on the symmetry axis is about 3 /inr, giving a contri¬ 
bution to the density variation uncertainty on the or¬ 
der of 0.0003 g.cm~ 3 . To this systematic uncertainty we 
add another incertitude due to the difference between 
the left and right parts of the Abel inversion. This gap 
varies from one image to another and it is on the order 
of 0.0005 g.cm~ 3 . 

We have measured many minimum densities in the 
vicinity of the cavitation voltage at several static pres¬ 
sures (0.15 bar, 0.65 bar and 1.26 bar) and at the same 
temperature 0.96 K as shown in Fig(T] In this figure, the 
error bars represent the quadratic sum of statistical and 
systematic uncertainties for each minimum density mea¬ 
surement. Within the error bars we find as expected 
that p C av is independent of P s t . Thus, we computed 
the mean of these measurements and their mean squared 
error to determine, respectively, the cavitation density 
of helium and its uncertainty. Our final result is that 
the cavitation density of superfluid helium-4 at 0.96 K 
Pcav = 0.1338 ± 0.0002 g.cm -3 . 

Note that the uncertainty given here (2 • 10 -4 g.cm -3 ) 
compared to the density variation Sp gives a relative un¬ 
certainty about 2%. Concerning the reliability of this 
measurement, it may be interesting to recall that a com¬ 


parison with a hydrophone [l] was made in water in 2010. 
It was found that the deviation between the two methods 
is less than 5%. 


VI. DISCUSSION 

Previous results @,H] about cavitation in liquid helium 
were given as cavitation pressures instead of cavitation 
densities. The equation of state (EOS) of liquid helium 
in its metastable state (density and pressure below the 
boiling curve values) is needed to convert the p cav to a 
corresponding P cav . Although such an equation of state 
has never been measured experimentally, some have been 
proposed. H. Maris has pointed out that, in the stable 
phase at T = 0.1 I\ , the sound velocity pressure depen¬ 
dence could be fit very well by the law c 3 = b(P — P c ) 
with c the sound velocity, P the pressure, P c the spin- 
odal pressure and b a constant [3 . He proposed that this 
relationship holds in the metastable state (negative pres¬ 
sure). Bauer et al. have performed Path-integral Monte 
Carlo simulations of liquid helium in the nretastable state 
at finite temperature and found the same dependence of 
sound velocity on pressure [3 • Dalfovo et al. have calcu¬ 
lated the EOS of metastable liquid helium at T = 0 K us¬ 
ing density-functional approach^, and Boronat et al. us¬ 
ing a quadratic diffusion Monte Carlo method to achieve 
a similar EOS[§j. The EOSs at 0 K agree within a few 
percent. Moreover, using the density-functional theory 
of Dalfovo et al., Maris and Edwards have shown that in 
the temperature range 0 < T < 1 K, the EOS is nearly 
independent of temperature [ill] . 

So in order to compare our cavitation density result 
to cavitation pressure results of other experiments, we 
use the well established EOS of metastable liquid helium 
at T = 0 K and assume it holds for T = 0.96 K. By 
doing so, our cavitation pressure is P cav (0.96K) = —5.1± 
0.1 bar. 

At temperatures ~1 K, in addition to the present ex¬ 
periment, there are to our knowledge only two exper¬ 
iments which studied the cavitation of liquid helium. 
Both also used focused acoustic wave. Xiong et al.\ 4] 
found the cavitation pressure at 1 K is ~ —3 bar. The 
incertitude mentioned in this paper is about ±10% and 
comes mostly from the difficulty of estimating the pres¬ 
sure at acoustic focus knowing the displacement of the 
emitter. Non-linear effects were not taken into account 
in their calculation. So this incertitude is likely to be 
underestimated. Caupin et al.fb] studied the dependence 
of cavitation voltage to the static pressure. They claim 
that this method enables them to set an upper limit 
for the actual cavitation pressure. Modelling a linear 
response of their emitter to voltage, they also give a 
lower limit for the cavitation pressure. Their result is 
—9.8 < P ca , v (0-9K) < —7.7 bar. According to the data 
points published in [15], the result at T ~1 K is almost 
the same. One can see that there are large discrepancies 
among these experiments. 
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We have tried to reproduce the experiment of F.Caupin 
et al. using their extrapolation method on P st [H] (see 
appendix). The upper limit of P cav we found is about 
—8 bar which agrees pretty well to the one of F.Caupin et 
al.. But the disagreement with our density measurement 
converted to pressure remains. 

Jezek et aZ.[l6) have calculated the cavitation pressure 
of liquid helium as a function of temperature, by using 
a density functional method and assuming the absence 
of defects (especially vortices). In order to compute the 
cavitation pressure, the volume v and the time r in which 
nucleation is likely to occur are needed. We take v = 
(A s /2) 3 and r = 0.1 /us is the 1/10 of the sound period. 
This gives vt ~ 10~ 13 cm 3 s. Using this vt value, Jezek 
et al. calculated P^ t f ek (0.96K) ~ —6.9 bar. This value 
is just between our result (—5.1 bar) and the central value 
(—8.8 bar) of reference 

Finally, we would like to point out that Maris has de¬ 
veloped a model of cavitation in the presence of quan¬ 
tized vortices in liquid helium G2- For a vortex density 
ranging from 10 4 to 10 12 cm -2 , he founds that —5.8 < 
^c" es (0.96IO < —5.1 bar. Although Maris can not 
estimate the error bar on this simulation, we note that 
our result does lie in this range. Besides, Pettersen et 
al. [18) have proposed that the vortex density in the 
high amplitude sound wave should be of the order of 
10 8 ~ 10 lo cm -2 . The presence of vortices might be a 
possible way to conciliate our experimental result with 
simulations. However, this would imply that, in the pres¬ 
ence of vortices, the P s t extrapolation method of refer¬ 
ence [1] does not give an upper limit of P ca v 


VII. CONCLUSION 


Using an interferometric set up, we have measured the 
cavitation density of liquid helium-4 at T = 0.96 K and 
the result is p cav = 0.1338 ± 0.0002 g/cm 3 . Trying to 
compare this result with existing calculations on the cav¬ 
itation pressure, we found that a model taking into ac¬ 
count the presence of vortices in the liquid can rather 
satisfactory explain our result. We plan to investigate 
the influence of vortices on cavitation in liquid helium 
in two ways. First by studying the dependence of p cav 
on temperature, a signature while crossing the lambda 
temperature should be seen. Second, we will probe the 
dependence on the density (in the metastable state) of 
the sound velocity and of the sound attenuation. This 
last part will be done using stimulated Brillouin scatter¬ 
ing. 
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Appendix: Cavitation pressure by extrapolating the 
static pressure 

In 2001, F.Caupin et al. implemented an extrapola¬ 
tion method [sj to investigate the behavior of helium in 
negative pressure. They imagined an environment with 
a “stable” negative static pressure. In this situation, the 
required driving voltage for achieving cavitation would 
be less than the one in null static pressure. Then the 
very negative static pressure corresponding to zero cav¬ 
itation voltage would be the cavitation pressure. That 
can be expressed as : 

Pcav = Pfocus = Pst + AP (PstVc ) (A.l) 

where P caV i Pst and Pf OC us are respectively the cavitation 
pressure, the static pressure of helium and the pressure 
at acoustic focus. A P ( p s tV ) is the variation of pressure 
induced by the sound wave and p s t is the static density. 
This equation holds when the driving voltage V reaches 
cavitation voltage V c . Assuming that the cavitation pres¬ 
sure is independent of P s t , we measure the different cav¬ 
itation voltages at different static pressures, and then 
extrapolate linearly at zero cavitation voltage (l9j| . Nu¬ 
merical simulations [uj) have shown that in absence of 
vortices the true curve is concave toward negative pres¬ 
sures. The linear extrapolation gives an upper limit of 
the true cavitation pressure. 
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FIG. 5. Reproduction of reference [§j experiment. P s t as 
function of p s tV c . The 3 data points(blac.k circles) correspond 
to static pressures of 0.15 bar, 0.65 bar, 1.26 bar. Cavitation 
voltage are the ones shown in Fig[2] The dotted line is the 
linear extrapolation of the date points. 

The Fig [5] shows our extrapolation corresponding to 
our measurements of ( P s t.,V c ) values (see Figj2|). The 
upper limit of cavitation pressure obtained in this way 
is —7.9 ± 0.3 bar. This is in agreement with F.Caupin et 
al.’s value. 
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